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Abstract: The highest energy cosmic ray event reported by the Auger Observatory has 
an energy of 148 EeV. It does not correlate with any nearby (z<0.024) object capable of 
originating such a high energy event. Intrigued by the fact that the highest energy event 
ever recorded (by the Fly's Eye collaboration) points to a faraway quasar with very high 
radio luminosity and large Faraday rotation measurement, we have searched for a similar 
source for the Auger event. We find that the Auger highest energy event points to a quasar 
with similar characteristics to the one correlated to the Fly's Eye event. We also find 
the same kind of correlation for one of the highest energy AGASA events. We conclude 
that so far these types of quasars are the best source candidates for both Auger and Fly's 
Eye highest energy events. We discuss a few exotic candidates that could reach us from 
gigaparsec distances. 
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1. Introduction 

Observatories such as the Pierre Auger and the High Resolution (HiRes) Fly's Eye have 
recently shed light on some of the open questions regarding the high energy end of the 
cosmic ray spectrum. However, much is yet to be revealed. Among the open questions 
are the origin of ultra high energy (above ~ 6 x 10 19 eV) cosmic rays (UHECR), the 
acceleration mechanism involved in their production, and their composition. 

As a first step in the search for the astrophysical sources of UHECR, the Auger col- 
laboration has observed an anisotropy of the UHECR arrival direction at 99% CL Q, via 
a correlation of cosmic ray arrival directions with nearby Active Galactic Nuclei (AGNs). 
The correlation with this anisotropic pattern in the sky does not imply that the AGN are 
necessarily the sources of UHECR, but rather that the cosmic rays tend to come from where 
the large scale structure of matter density is distributed. It suggests that the deflection of 
cosmic ray trajectories by galactic and extragalactic magnetic fields may be small enough 
that cosmic ray astronomy becomes possible. To add to the puzzle, a similar correlation 
with nearby AGN has not been seen in the HiRes data ||. 

Both experiments have observed (3|, U the expected Greisen-Zatsepin-Kuzmin (GZK) 
suppression Q in the cosmic ray spectrum due to resonant scattering on the cosmic mi- 
crowave background photons. Furthermore, the Auger correlation with nearby AGN is 
compatible with the expectation of the resulting GZK particle horizon. 
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Here we focus our research on the highest energy events, those with energy above 
10 20 eV. Naively, one would expect that the highest energy UHECRs would be the least 
deflected by magnetic fields, and hence would be the events most likely to point back to 
their astrophysical sources. Paradoxically, there are no known plausible sources within 
the GZK sphere for the highest energy events, leading some to hypothesize that UHECR 
are produced by nearby, energetic, transient sources such as gamma ray bursts B, |7| or 
AGN flares || |9[. These sources would not necessarily have large photon fluxes observed 
concurrently with a UHE hadronic cosmic ray flux. 

These events might also be heavy nuclei being deflected by the Galactic magnetic field. 
As it is easier to accelerate heavy nuclei than protons, one can analyze the possibility of 
heavy nuclei acceleration. 

Another exciting possibility is that new, exotic physics may be responsible for the 
unimpeded propagation of UHECRs over cosmological distances. Although the presence of 
the GZK feature is compatible with a standard particle composition of UHECR, it does not 
rule out the possibility of a small fractional flux of an exotic component. Such an exotic 
flux may indeed be able to penetrate through the opaque wall of background photons and 
travel far larger distances without attenuation. 

Inspired by a thorough search [1C] for the source of the highest energy event so far 
detected (by the Fly's Eye collaboration fll|]), we search for a possible source of the Auger 
highest energy event. This event does not correlate with any nearby AGN from the Veron- 
Cetty and Veron catalog [^] used in the Auger correlation studies nor with any nearby 
source capable of producing such an energetic event. This lack of correlation is particularly 
puzzling because the event points back to an apparently well-covered region of the sky 
between the Virgo and Centaurus superclusters, where the VC catalog would be expected 
to be complete for the highest luminosity, nearby AGN. Instead, we find that the most 
energetic event so far detected by the Auger collaboration, points to a faraway quasar, 
well beyond the GZK radius. This quasar has the same characteristics as the one which 
correlates with the Fly's Eye event. 

We also analyze the arrival direction of AGASA's 11 events above 10 20 eV [14|. It 
should be noted however that AGASA's energy scale is systematically 30% higher than 
the HiRes/Fly's Eye energy scale [13] while Auger's energy scale is around 15% lower. 
Unlike the Auger and HiRes measured UHECR spectra, the AGASA spectrum at high 
energies exhibits no GZK suppression; this may be indicative of poorer energy resolution 
at these energies. The HiRes collaboration has not yet published their event list containing 
the energy and direction information together. For this reason, it is not possible for us 
to analyze their results as we did for Auger and AGASA. However, it is intriguing that 
HiRes has already reported a potential correlation of their stereo events with distant BL 
Lac objects, well beyond the GZK horizon [18, |17]]. This correlation cannot be explained 
by conventional physics models. 

In the next section we will review the search for the source of the Fly's Eye highest 
energy event. In section |] we describe our search for source of the Auger highest energy 
event. Then follows the analysis for AGASA's events. In section |6] we discuss exotic models 
which may allow cosmic rays to reach us from such faraway sources. 
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2. The highest cosmic ray event ever detected 



The highest energy cosmic ray event was detected by the Fly's Eye collaboration [O], which 
determined the energy to 3.2 x 10 20 e V. A thorough search for its source was done by Elbert 
and Sommers (ES) p0|| , who looked for strong radio sources with strong magnetic fields 
and large fluxes of energetic particles in the direction of the Fly's Eye event. These sources 
were targeted since in order to accelerate particles to ultra high energies, a high magnetic 
field moving with high velocity is required ]19|| . Powerful radio galaxies are therefore good 
UHECR source candidates. Their search was limited to a region within an angular distance 
of 10° from a two sigma error box centered in the event arrival direction. 

From the four sources that fell into this region, the quasar 3C 147 was the only one 
within one sigma from the event arrival direction. It is also the only one which meets the 
necessary requirements for particle acceleration. It is listed in a catalog of 173 powerful 
radio luminosity sources ^0|. Its total radio luminosity is about 10 45 ergs s _1 being among 
the 28 brightest ones in this catalog. It is also included in a list of 96 plausible candidates 
for having very large Faraday rotations (FR) |21|, having one of the top four largest FR 



22]. The redshift of this source is however 0.545, which puts it at a distance of 2 Gpc away 



from us for current cosmological parameters. 



3. Procedure for Source Correlation 

In this section we describe our procedure for correlating sources with the arrival direction 
of the high energy events. We follow the same procedure adopted by ES in looking for high 
luminosity sources in the direction of the Fly's Eye event. In order to correlate the HEA 
to a faraway source, we assume the cosmic ray is electrically neutral. For this reason, the 
angular cut around the event arrival direction is defined as 1.5 x a exp , where a exp is the 
experimental angular resolution. 

As for the Fly's Eye event, we are not looking for any object in the direction of the 
detected event, but for sources with the capability of accelerating particles to 10 20 eV. 
It has been shown p3"fl that a solid lower bound on the source luminosity in order to 
achieve such an acceleration, is of 10 45 Z -2 erg/s, where Z is the charge of the cosmic ray 
nucleus. Our goal is to find sources with 10 45 erg/s minimum luminosity in the direction 
of the detected event. We should point out however that although this bound is valid for 
many acceleration models including the diffusive or Fermi mechanisms, it does not hold for 
mechanisms like inductive acceleration that might occur near central black holes in AGNs. 
Also the requirement is on the magnetic luminosity but as in [^] we use it as a bolometric 
luminosity limit. 

As a first step we select high luminosity sources according to a minimum flux density. 



We use the high luminosity criterion established by Wall &: Peacock [24], which requires 
a minimum flux density of 2 Jy at a frequency of 2.7 GHz. This allows us to produce a 
list of astrophysical objects, which although broad is a first selection step. Once this list 
is produced we verify which of these objects have the minimum required luminosity. 
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In order to account for measurement uncertainties, we require a minimum source lu- 
minosity of L = 10 44 erg/s. One can roughly translate this to a minimum flux density -Smm 
at a certain frequency v in Hz [|lOf |: 

o ^ LHj 5xlQ9 
Attuc^z^ UZ z 

where c/Hq is the Hubble length, z is the source's redshift, and S^nm is given in Jansky. We 
use the luminous parameter Rl = S/S m in defined by ES, where S is the source flux density 
in Jansky and require Rl > 1- As an example, the nearby Centaurus A has Rl = 10 -3 
and M87, 2 x 1CT 3 . As we will see the quasar 3C 147 found in the direction of the Fly's 
Eye event has Rl = 6. This will reduce the number of objects selected and allow us to 
require a minimum angular distance from the event arrival direction and then integrate 
the observed flux density to determine the luminosity. The total radio luminosity can be 
determined as [^0| 

L = 4 X 2 (l + z) 2 S bol , obs h- 2 (3.2) 

where x ls the comoving distance and Sboi,obs is the observed radio spectrum integrated 
from 10 MHz to 100 GHz. 

Both the first list of high luminosity objects as the flux density at various frequencies 
can be obtained from the NASA/IPAC Extragalactic Database (NED) catalogue p5fl . 
For each telescope we look for sources within their field of view. Extra requirements are 
that the sources are away from the galactic plane (|6| > 10°) where astrophysical catalogs 
are incomplete, and located beyond the GZK radius, at a redshift larger than 0.024. The 
GZK radius is the maximum distance from which a standard particle can travel through the 
cosmic microwave background radiation (CMB) and reach the Earth at ultra high energies. 
For a proton this distance is about 100 Mpc. The latter requirement comes from the fact 
that searches within the GZK cutoff were already done. The result for a nearby search for 
the Fly's Eye event source is described in section and we will describe the Auger search 
in section ||. It is worth noting that no high luminosity nearby source was found in the 
direction of these events. 

We now describe how we determine the chance probability that a specific correlation 
comes from an isotropic distribution. The probability p that the arrival direction of a 
particular event, which comes from an isotropic distribution of sources in the sky, points 
by chance to one specific source in the candidates list is given by: 

p= e(S) cos (6) d5 da ji (3-3) 

J A 

where the integration goes over an area A and e(d) is the experimental exposure function 
which depends on the declination 5 and the zenith angle 9 and is given by 

f 2 

e(8) cx / 7r cos OdoiR (3-4) 
Jo 

where or represents the right ascension and 

cos 9 = sin 5 sin / + cos S cos I cos or (3-5) 
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where I represents the observational site latitude. The exposure function is normalized to 
unity over the experimental field of view of the observatory. The area A is the area in 
the sky determined by circular regions of radius (5 around each candidate source. For the 
Auger and Fly's Eye telescopes, there is no exposure for zenith angles above 60° and the 
integrand of equation [O] is zero in this region. For AGASA there is no exposure for zenith 
angles above 45°. 

We determine p through a Monte Carlo simulation, where we generate a large number 
of sources isotropically distributed over A, where they follow the pattern determined by the 
exposure function e(5). We then compute the angular distance between these simulated 
events and the sources that belong to our catalogue, p is determined considering that the 
event arrival direction matches the candidate source direction within an angular distance 
of 1.5a exp . 

3.1 Reanalyzing the Highest Energy Event 

Here we reanalize the the Fly's Eye highest energy event, looking for bright sources in its 
direction. We follow each of our procedure steps described above. We aim at the same 
time to check our method and if no other source - that could have been discovered after 
the work of ES - also matches the arrival direction of the Fly's Eye event. 

One important point in this analysis is that the Fly's Eye telescope exposure is a 
complicated function of energy and zenith angle. For this reason our analysis determines 
the Fly's Eye exposure in a very crude way. The Fly's Eye angular resolution is asymmetric 
in right ascension and declination, being the later much worse with a value of 6 degrees. To 
simplify our search we look for sources within 9 degrees from the event direction. This will 
of course make the chance probability of having a correlation from an isotropic distribution 
very high. 

We produce our first source list by requiring a minimum flux density of 2 Jy at a 
frequency of 2.7 GHz. This list contains 148 sources. The second step is to require a 
minimum flux density as described by Rl, which reduces this list to 54 objects. 

Our result is compatible with the one found by ES. We find that the only high lumi- 
nosity object that correlates with the Fly's Eye event direction is the Quasar 3C 147. The 
caracteristics of this quasar is described in section || and also in the work of ES. The chance 
probability for this correlation is large, being p = 0.37 if we consider all 54 sources in our 
final catalalogue. ES quote that the number of sources from their catalogue expected to 



fall in their 2a error box centered in the event direction is 0.24. We integrate equation 3.2 
and determine 3C 147 luminosity to be 10 erg s _1 , the same as found by ES. 

This shows that our procedure is compatible with the one done by ES and we now use 
it to analyze both Auger and AGASA's events. 



4. The Auger Highest Energy Event 

From the 27 high energy events published by the Auger collaboration Q, only one is above 
10 20 eV, having an energy of 148 EeV. Events of these energies are expected to have its 
origin outside of our Galaxy fl~9|| . As ordinary galaxies do not meet the requirements for 
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accelerating particles to such high energies, we searched for high luminosity radio sources 
in the direction of the highest energy Auger event (HEA). We proceed as described in 
section ||. 

We produce our first source selection list requiring a minimum flux density of 2 Jy 
at a frequency of 2.7 GHz in the "advanced all-sky survey" of the NED. The equatorial 
coordinates for the HEA are right ascension or = 192.7° and declination 5 = —21°. The 
Auger exposure Jjj] for events with zenith angle below 60° is roughly constant in right 
ascension and is a function of declinations, being nonzero for 5 < 24.8°. 

The NED search results in 184 strong radio luminosity objects. We further reduce this 
list by requiring £ m i n as defined in Equation [O] which is equivalent to Rl > 1. From the 
184 NED sources, only 67 meet this minimum luminosity requirement, among which 57 are 
quasars and 10 are galaxies. As described in section |3| we look for sources with a redshift z 
beyond 0.024. As the Auger optimized search adopted z < 0.018, we also look for sources 
with a redshift between 0.018 and 0.024. This search is described a bit further. 

As for these high energies the Auger angular resolution is about 1° Q, we require that 
the angular distance between the source and the event is at maximum 1.5°. From our 67 
objects catalog, only one fell in this angular region. It is the quasar (QSO) PKS1245-19, 
with an angular distance of 1.2° from the HEA. The chance probability for this to happen 
from an isotropic distribution is 0.013 1 . 

The QSO PKS1245-19 is also included in a sample of the most powerful radio sources 



from the Molonglo catalog [27]. Molonglo is a 408-MHz survey looking for sources with 
flux densities above 0.7 Jy. This catalog covers 7.85 sr of the sky, over declinations 
+18°. 5 > 5 < — 85°. and \b\ > 3°. No other powerful radio sources are found in this 
catalog in the direction of the HEA 2 . PKS1245-19 is also among 96 sources selected as 
good candidates for having high Faraday Rotation [21]. It has also been listed as an 
energetic gamma-ray source detected by the Egret telescope [pO] . 

This source candidate is therefore very much like the one found in the Fly's Eye event 
direction (QSO 3C147). An integration of the PKS1245-19 radio spectrum (equation |3.2| ) 
gives a radio luminosity of ~ 3 x 10 46 erg/s. As both QSOs are also listed as having high 
Faraday Rotation measurement, they are excellent candidates for accelerating particles to 
ultra high energies. Just as the case with 3C147, PKS1245-19 is very far away, at a redshift 
of 1.275 (3.8 Gpc for current cosmological constants). 



1 Although efforts were made to fix the analysis before looking at the data, mistakes were made and the 
data was processed through two different sets of cuts on b and the maximum angular distance between 
the source and the cosmic ray arrival direction before the analysis was finalized. So at most, the reported 
probability would be multiplied by a factor of 3 to account for previous trials. 

2 If we modify our search to look for objects within z < 0.024, we find 19 sources which fulfill our initial 
requirement of at least 2 Jy at a frequency of 2.7 GHz. From these only one has an angular distance a of 
less than 15° from the HEA. It is the galaxy pair VV201 [0 at z = 0.015 or a distance of about 61 Mpc, 
and q = 8.5° from the event direction. If a 1 Mpc magnetic field coherence length is assumed, an rather 
large extragalactic magnetic field of ~ 4 nG is necessary to bend a proton to the event direction [^9j. The 
radio luminosity of this pair of interacting galaxies is ~ 5 x 10 41 erg/s, much lower than the 10 45 erg/s 
required for UHECR sources. This source does not fulfill the S m in requirement, Eq. p.l|. 
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5. AGASA High Energy Events 



The Akeno Giant Air Shower Array (AGASA) has detected 11 events with energies above 
10 20 eV. Although there is a ~ 30% discrepancy among their energy measurement and 



HiRes' [13], we take their results as published |L4|]. Given these uncertainties we search 
for far away sources correlating with the AGASA events for completeness. There are other 
works searching for sources of AGASA events [15, [l6|] which find no significant correlation. 

Following our prescription we search for very high luminosity (above 10 45 erg/s) objects 
in the direction of any of these 11 events. The only difference is that due to the poorer 
AGASA angular resolution, we look for a high luminosity source within a 3° angular 
distance from each event direction. This will of course increase the chance probability of 
having an incorrect source in the direction of an event, if sources positions are sampled 
from an isotropic distribution. The AGASA latitude is 35°47' and its field of view covers 
approximately between —10° < 5 < 80° in Equatorial coordinates. Our final catalogue of 
high luminosity sources in the AGASA field of view consists of 54 objects. From these the 
QSO 3C380 is at an angular distance of 2.5° from AGASA's 1.05 x 10 20 eV event. This 
event has a = 298.5° and 5 = 18.7°. The chance probability for a 3° correlation from an 
isotropic distribution for one isolated event correlation is of 0.048, and accounting for the 
fact that one event among 11 correlates makes this probability equal to 0.42. 

Estimating the radio luminosity of this quasar from Equation |3.2| , one gets ~ 10 46 erg/s. 
It is also listed in the same large Faraday Rotation list [^] as the above QSOs. Just as the 
previous candidate sources, QSO 3C380 is an excellent source for acceleration of UHECR. 
It is also faraway, at a redshift of 0.692 (2.4 Gpc for current cosmological parameters). 

Since the chance probability for source correlation with AGASA's events directions is 
very large, we might conclude that there is no significant correlation with high luminosity 
objects. However as discussed in the introduction there are uncertainties on the energies 
of these events and they might not fulfill our requirement of 10 20 eV minimum energy. 



6. Exotic candidates 

As described in the previous sections, QSOs PKS1245-19, 3C147 and 3C380 whose positions 
are well-correlated with the directions of Auger's, Fly's Eye's highest energy events, and to 
one of AGASA's high energy events, are very good UHECR source candidates. The only 
question then relates to their large distance from the Earth. 

One possibility, that involves new physics but not necessarily new particles, is that the 



propagation distance of UHECR can be increased if Lorentz symmetry is violated [31]. In 
this work we analyze possible new particles that could be responsible for the high energy 
events. 

Below we describe two candidates that could account for this huge distance. One 
important remark is that although the GZK feature has been observed, high energy exotic 
events can certainly compose a small fraction of the UHECR spectrum. They could maybe 
even account for most of the composition at the highest end of the spectrum above 10 20 eV. 
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Two possible candidates would be axions-like particles [32] and exotic massive hadrons [33 



34]. Photons would convert into axions at the source, due to their interaction with the 
source's magnetic field. They would then propagate through intergalactic space without 
scattering on the CMB or other extragalactic background photons, and convert back to 
photons in the Milky Way. Their air shower would thus be initiated by a photon. How- 
ever, as we will describe below, the photon-axion conversion probability is low at such 
high energies. Alternatively, exotic massive hadrons would be the primary particle in the 
atmospheric shower development |54|. These showers can be distinguished from showers 
initiated by standard particles [^5|. However, as we will describe later, it is not clear how 
these particles can be produced without being constrained by the diffuse photon flux, since 
a large number of photons should be produced simultaneously. 

6.1 Axion-Like Particles 



Axion-like particles (ALP) []32| would be excellent candidates to propagate through cosmo- 
logical distances without scattering. They might also explain the possibly anomalous trans- 



parency of the universe to TeV gamma rays observed in air Cherenkov telescopes [37, 38]. In 
models of photon-axion mixing, an UHE gamma is produced as a byproduct of acceleration 
of hadronic cosmic rays in the astrophysical source. This gamma coherently interacts with 
background magnetic fields and is converted into a UHE axion which can then propagate 
cosmological distances without scattering. Upon entering the magnetic field of the Milky 
Way, the axion reconverts into a UHE gamma which induces a detectable cosmic ray air 
shower. Due to low statistics at the highest energies, experimental constraints on a possible 
gamma-induced fraction of observed air showers are still rather weak, with upper bounds 
ranging from the 20% — 60% level for UHECR of energies greater than a few 10 19 eV (see 
for example [fDjI ). It has also been shown that the longitudinal profile of the most energetic 
Fly's Eye event is compatible with that expected of a photon-initiated air shower at the 



1.5a level p2 [. Morever, due to the poorer composition sensitivity in UHECR detected by 



surface detectors |39|, |40|], it is not possible to unambigously distinguish between hadronic 
or photon-initated showers for either the Auger event or the AGASA events. In particular, 
there is no composition-sensitive muon flux information on the correlated AGASA event. 
So it is possible for all three of the correlated events discussed here to be photon-initiated. 
The interaction term between the axion field a and photons is 

C- mt = ~gaFF = gE ■ B (6.1) 

where g is the axion-photon coupling, F represents the photon field, and E and B represent 
the electric and magnetic fields respectively. In the presence of a background magnetic field 
of strength B, this interaction term becomes an off-diagonal term in the effective flavor 
mixing matrix, coupling axions to the photon polarization component aligned with the 
magnetic field direction. At sufficiently high photon energies u > (m„ — m^)/(2gB), the 
photon-axion system is in a regime of maximal mixing, and the probability of conversion 
from photons to axions and vice- versa is P ~ {gBL) 2 /A. This probability becomes of order 
unity in regions where BL ~ 1/g. 
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The UHE photons are expected to be produced as secondary particles in the accleration 
of UHE protons. In order to accelerate cosmic ray protons to -Bproton = 10 20 eV energies, 
the sources must have magnetic field regions which satisfy BqL = E pTOtOQ . This coincides 
with the results of galactic magnetic field models |3(| which arrive to a comparable magnetic 
baseline for the poloidal magnetic field, BqL ~ (6^xG) • (4 kpc) = 3 x 10 19 eV. So an axion 
model with g ~ 10~ n GeV -1 would predict efficient conversions via oscillation both at the 
source, and in the galaxy [37, |3£|. Alternatively conversions could occur in propagation 



over sufficiently large coherent domains of the weaker extragalactic magnetic fields [ p3j| . 

However when considering photon-axion oscillations at ultra high energies, effects from 
magnetic birefringence might dominate implying in an insignificant oscillation probability 



1 44]. At energies E » TeV, magnetic birefringence due to quantum electrodynamics 
alters the photon dispersion relationship and induces a large diagonal component into the 
effective mixing matrix. As a result the effective photon-axion mixing angle becomes too 
small to produce efficient photon-axion conversion. A previous paper proposing photon- 
axion mixing as the source of the correlations between HiRes cosmic ray events and distant 
BL Lac objects [17] did not account for the QED birefringence effect 45], so here we will 
investigate the effect in detail. 

The formulation of the birefringence of the vacuum, was first quoted by Adler [46], 



following numerical evaluations by Toll [47]. In the presence of a constant background mag- 



netic field B, the indices of refraction for photon polarizations perpendicular and parallel 
to the B are given by: 

a ( BsinO^ 2 



where a is the fine-structure function, x parameterizes the center of mass energy as 

u B sin 9 



2m B r 



(6.3) 



9 is the angle between the photon propagation and the direction of the background field 
B, and the critical magnetic field is 



B CI = m z Je = 4.41 x 10 id G. (6.4) 

The function N±u(x) is determined by Toll [47], by numerically evaluating the dispersion, 
averaged over a series of absorption resonances in which the external magnetic field absorbs 
the excess energy-momentum in the process 7 — > e + e~ . The numbers given by Adler, and 
subsequently widely quoted, including in [48] are 3 



N± =8/45, (6.5) 
Nu = 14/45. (6.6) 



3 Note that we use here the conventional definition of photon polarization using the electric field vector 
rather than the opposite definition used by both Toll and Adler. 
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Most subsequent astro-physical photon-axion mixing papers [49, 50, 51, 52 1 primarily 
concern themselves with mixing at low energies, and therefore do not explicitly quote 
Adler's admonition that these values only represent a particular asymptotic limit of Toll's 
calculations when B <C B c and ui < 2m e . This is the limit when the available energy in the 
interaction between the photon and magnetic field is far below the first pair production 
resonance, a regime of normal dispersion. In Toll's thesis, it is illustrated that as x increases 
through the series of pair production resonances, whose effects can be averaged over, the 
photon enters a regime of anomalous dispersion where the index of refraction decreases 
through n = 1 before becoming slightly negative. In the limit B <C B cr , the numerical 
predictions of the asymptotic behavior at high x far beyond the resonances are 



N ± 



-0.4372/x 4 / 3 , 
-0.6558/x 4 / 3 . 



(6.7) 
(6.8) 



This approximation is valid when x > 500 or B m Q u\2 > 10 13 , where magnetic field is 
measured in milliGauss, and photon energy in units of TeV. 
The axion-photon mixing matrix is given by pg|] 
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(6.9) 



Where A±_ and An are orthogonal components of the photon field. The elements of the 



dispersion matrix are defined and estimated [44] as 
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(6.10) 
(6.11) 

(6.12) 

(6.13) 



where m^y = m a /fieV, gn = g x 10 11 GeV, is the plasma frequency for an electron 
density n e . As we have now seen, the last equation holds only for low magnetic field and 
low energy photons such that x <C 10 _1 or B m Q ■ u\2 10 10 . 

Since the sign of the QED birefringence term Ajj^ ED is negative, indicating a regime 
of normal dispersion, and A p i is negligible, there is no possibility to achieve the resonance 
condition A^ ED = A a where the mixing angle 6 = (1/2) arctan(2AB/(AQED + A p i — A a )) 

becomes maximal. The best that can be hoped for is that both A^ ED ,A a Ab so the 
configuration is sufficiently close to resonance that a large mixing angle can still be achieved. 
Suppose that the initial conversion from photons to axions occurs in the galaxy hosting 
the astrophysical photon source, and the regeneration of photons occurs in the Milky Way 
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galaxy. Then, the highest energy photons that can penetrate the wall of extragalactic 
background light can be estimated using g\\ = 10, B m Q = 10~ 4 , and L ~ 30 kpc for both 
source and destination galaxy so that the mixing probability is of order unity in each. 
In this case, the maximum photon energy for which this process can occur efficiently is 
10 17 eV, for which A a ~ 10" 25 cm" 1 , A^ ED ~ A B ~ 10~ 24 cm" 1 . 

The only remaining possibility to obtain efficient mixing at even higher energies is 
to enter the regime of anomalous QED dispersion, in which the indices of refraction are 
parameterized by Eqs. ^3~J. In this case, A^ ED > 0, and so an exact resonance is possible. 
To achieve this condition, the center of mass energy must be well above the energy required 



to produce on-shell pairs, B m Q ■ oj\2 > 10 | 47fl . For photon energies of 10 eV, the 



magnetic field required is then 10 4 G. Such high fields may indeed be present in the high 
energy photon sources, and induce resonant conversion of photons into axions, or even 
efficient conversion by adiabatic level-crossing. However, there are no such fields in the 
Milky Way, and the resulting ultra-high-energy axions cannot be efficiently reconverted 
into detectable photons. 

6.2 Exotic Massive Hadrons 

Before the high energy end of the cosmic ray spectrum was established, a class of exotic 



massive hadrons was proposed [33] to account for the events beyond the GZK cutoff. This 
generic class was coined as "uhecron", and is composed of colored, strongly interacting 
massive particles. 

Uhecrons are modeled ]34j as a heavy single constituted core surrounded by light 
hadronic degrees of freedom. The core concentrates the momentum of the particle while 
the light degrees of freedom (which can be gluons and/or light quarks) is responsible for the 
interaction. The strong interaction is necessary in order to produce a hadronic shower in the 
atmosphere [|34| . If the uhecron is too heavy, the momentum carried by the light consituent 
will be too small since most of the energy is carried by its core and no interaction in the air 
will occur. It was shown |34j] that in order to produce a hadronic shower in the atmosphere, 
an uhecron cannot be heavier than 50 GeV. It is important to note that particles such as 
the proposed S° Q cannot be classified as an uhecron, since its momentum is shared 
almost equally by all its constituents. 

Uhecrons can avoid the GZK cutoff since it loses much less energy than protons while 
traveling through the CMB. In short, their heavy mass shifts the threshold for photopro- 
duction of pions to much higher energies. They can propagate distances of gigaparsecs 
through the CMB before losing most of its energy |33| . Uhecrons are also taken as electri- 
cally neutral particles in order to avoid deflections by magnetic fields. 

Simulations of uhecron induced air showers |M| show that uhecrons with masses up 
to 50 GeV are compatible with UHECR events. The shower generated by the interaction 
of an uhecron with the atmosphere, can resemble the longitudinal profile of the highest 



energy cosmic ray detected by Fly's Eye [ljj . 



One of the best candidates for an uhecron is the heavy gluino lightest supersymmet- 
ric particle (LSP), proposed by Stuart Raby |53|]. It fits all the uhecron requirements. 
Experimental limits set the heavy gluino LSP mass between 25 and 35 GeV [p^j. 
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As uhecrons can account for the large distances between faraway QSOs and the Earth, 
the heavy gluino LSP is a candidate for the highest energy cosmic rays analyzed here. 

It has also been showed that it is possible to distinguish uhecron from proton or heavy 
nuclei induced showers [pqj . Uhecrons with masses in the 25 to 35 GeV window can be 
discriminated from the bulk of the cosmic ray spectrum, which might be composed by 
protons or iron. An enhancement of this discrimination can be achieved by correlating the 
event direction with the far away source. 

However, the acceleration mechanism to produce uhecrons at these energies is not well 
established. It has been pictured generically [33] as originating from proton collisions with 
hadronic matter in the source. This requires protons being accelerated to energies slightly 
above 10 20 eV. It was pointed out |33l] that even if uhecrons are originated in proton 
collisions, it is hard to explain the lack of lower energies uhecrons. This remark does 
not hold here, since uhecrons would constitute a small fraction of the UHECR spectrum, 
and the lower energy ones would be hidden by the larger fraction of standard particles. 
Information from the Large Hadron Collider on the heavy gluino LSP will contribute to 
the possibility of uhecrons accounting for the long distances between QSOs and Earth. 

The main potential problem with uhecron production from higher energy proton col- 
lisions is that since the uhecron is modeled as a neutral particle, photons will also be 
produced in abundance. If all similar sources as the ones found both in the Auger and 
Fly's Eye event direction, produce uhecrons, a significant flux of photons should be found. 
As the Fermi measurement on the diffuse extragalactic gamma flux [55] is even lower than 
the one from Egret, it seems that most top-down scenarios are excluded as significant 



contribution for the UHECR spectrum [56 



7. Conclusions 

We have described our search for the source of the highest energy cosmic ray detected 
by the Pierre Auger Observatory. We show that given requirements on the acceleration 
capabilities of the source, one faraway quasar is located at 1.2° from the event arrival 
direction. The chance probability for this correlation from an isotropic distribution is of 
0.013. 

In this search we look for sources with radio luminosities above 10 45 erg/s, which are 
required in order to be able to accelerate particles to energies of 10 20 eV. This requirement 
is valid for many acceleration models, as the diffusive or Fermi mechanisms, but not for 
some mechanisms as the inductive acceleration close to central black holes in AGNs. Also 
although the bound is set for magnetic luminosities we use it as a bound for bolometric 
luminosities as in ]23[| . The quasar PKS1245-19 found in the direction of the Auger event, 
has a luminosity of 4 x 10 46 erg/s. It is included in a catalogue of the most powerful radio 



sources [27 and among 96 sources selected to be good candidates for having large Faraday 



Rotation [21|. It is also listed as an energetic gamma-ray source by the Egret telescope 



1 3C ] . PKS1245-19 is therefore very much like quasar 3C147 found |Kj at one sigma from 



the direction of the highest energy cosmic ray ever detected [11] 
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As no obvious astrophysical source within the GZK cutoff was found in the direction 
of either the Auger or the Fly's Eye event, the faraway quasars are so far the best observed 
candidate sources for these high energy cosmic rays. In this case the only explanation for 
these particles being able to reach us, is that they are composed of non-standard particles, 
at least during their transit, or a non standard propagation mechanism. 

Two plausible candidates are axions and exotic massive hadrons. As was shown in this 
work, the probability for photon-axion conversion at these high energies is probably too 
small for this effect to produce a significant contribution to the UHECR flux. Uhecrons 
could account for these events, but there are potential problems related to their origin. 
Accelerating protons to even higher energies above 10 20 eV is required in order to originate 
uhecrons at the required energy. This mechanism will simultaneously generate high energy 
photons, which flux is limited by current gamma telescopes. One should however bear in 
mind that the acceleration mechanism for standard particles is not well understood either. 
Furthermore, it could be that the energies of these exotic uhecron-induced showers are 
systematically overestimated since the energy scale of the various detectors is calibrated to 
the average, presumably hadronic, cosmic ray. If either the shower profile or the invisible 
energy carried by neutral shower particles differs from that expected from typical hadronic 
showers, then the energy measurement can be biased |)7|] for a small fraction of the cosmic 
ray flux composed of exotic particles. 

As we are analysing events at energies well beyond that produced in Earth based 
accelerators, there could still be some unknown exotic phenomena involved in the cosmic ray 
production and propagation. It has been pointed out that violation of Lorentz symmetry 
could greatly increase the propagation distance of UHECR by changing the kinematics 



of the scattering on background photons [31]. Although we cannot find non problematic 



candidates that could account for these events, it is clear that there exist faraway quasars 
in the direction of these ultra high energy events and that only new physics can account for 
a particle propagating from their distance to us. The alternative of course, is that UHECR 
sources are nearby, but transient and not readily apparent in current observations of the 
sky. Or that they are heavier nuclei which do not point to their sources. 
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